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Abstract 

We evaluate the neutrino fluxes to be expected from neutrahno LSP annihilations 
inside the Sun, within the minimal super symmetric extension of the Standard Model 
with super symmetry-breaking scalar and gaugino masses constrained to be univer- 
sal at the GUT scale (the CMSSM). We find that there are large regions of typical 
CMSSM (?Tii/2;"io) planes where the LSP density inside the Sun is not in equilib- 
rium, so that the annihilation rate may be far below the capture rate. We show 
that neutrino fluxes are dependent on the solar model at the 20% level, and adopt 
the AGSS09 model of Serenelli et al. for our detailed studies. We flnd that there 
are large regions of the CMSSM (?tii/2)"^o) planes where the capture rate is not 
dominated by spin-dependent LSP-proton scattering, e.g., at large along the 

CMSSM coannihilation strip. We calculate neutrino fluxes above various threshold 
energies for points along the coannihilation/rapid-annihilation and focus-point strips 
where the CMSSM yields the correct cosmological relic density for tan/3 = 10 and 
55 for ^ > 0, exploring their sensitivities to uncertainties in the spin-dependent and 
-independent scattering matrix elements. We also present detailed neutrino spectra 
for four benchmark models that illustrate generic possibilities within the CMSSM. 
Scanning the cosmologically-favored parts of the parameter space of the CMSSM, 
we find that the IceCube/DeepCore detector can probe at best only parts of this 
parameter space, notably the focus-point region and possibly also at the low-mass 
tip of the coannihilation strip. 
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I. INTRODUCTION 



The most convincing way to verify directly the existence of astrophysical cold dark mat- 
ter particles would be through their scattering on nuclei in low-background underground 
experiments [ij. Complementing this search for dark matter scattering in the laboratory 
via spin- independent and/or -dependent scattering, the next most direct way to search to 
confirm the nature of any such dark matter particles would be to observe products of their 
annihilations in an astrophysical context. Although less direct, observations of such anni- 
hilation products would provide valuable insight into the annihilation rates into important 
channels, which would provide more information on the dynamics of the dark matter par- 
ticles, beyond their scattering cross sections. Observations of astrophysical annihilations of 
dark matter particles would be particularly interesting because cosmological annihilations 
earlier in the history of the Universe controlled the primordial relic abundance. 

Several environments for annihilations of astrophysical dark matter particles are of inter- 
est, including the galactic halo j2|, the galactic centre jsj, the Sun the Earth [sl, dwarf 
galaxies jsj and galaxy clusters our focus here is on the flux of neutralinos from dark 
matter annihilations inside the Sun. In this case, the dominant astrophysical uncertainties 
are the total line density of cold dark matter through which the Sun has passed throughout 
its history, and the mass fractions of different elements in the Sun, particularly for heavy 
nuclei (metals). In this paper we do not discuss the former, but we compare the predictions 
of models with different element compositions and discuss the importance of this model 
dependence for the observability of a signal. 

There have been many studies of the prospective neutrino fluxes from solar dark matter 
annihilations in generic versions of the Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM), and also in specific versions such as the CMSSM [8|, in which the 
supersymmetry-breaking scalar and gaugino masses are each constrained to be universal at 
some input GUT scale, and the lightest neutralino, Xi is assumed to be the lightest super- 
symmetric particle (LSP) j9-ll|. In this paper we stud y y annihilations in the CMSSM, 
focusing in particular on the strips of parameter space [lO] where the relic x density lies 
within the range of cold dark matter density indicated by WMAP and other cosmological 
measurements. We also present more complete results for a few specific benchmark CMSSM 



scenarios 



13| . which are representative of the range of possibilities within the CMSSM. 
The first step in calculating the flux of neutrinos from LSP annihilations inside the Sun is 
to calculate the LSP capture rate. This is controlled by the cross sections for LSP scattering 
on protons, which is mainly spin-dependent in the CMSSM, and on heavier nuclei, which 
is mainly spin- independent. Early work on capture often assumed that spin-dependent 
scattering would dominate, and neglected the spin-independent contribution from heavier 
nuclei, in particular. Secondly, one must check whether it is correct to assume that the LSP 
capture and annihilation processes in the Sun have reached equilibrium. 

As we show later in this paper, neither of these assumptions is valid in general in the 
CMSSM: spin-independent scattering is also important, particularly at large values of mi/2, 
and the annihilation rate is significantly smaller than the capture rate in large areas of 
the (mi/2,mo) plane. Both these effects are particularly important at large mi/2 along 
the coannihilation WMAP strip. We also discuss the differences in the annihilation rates 
calculated using different solar models. 

The spin- dependent and -independent dark matter scattering rates both have significant 



uncertainties that affect the capture rate and hence also the annihilation rate |14j-|l8j. The 
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matrix element for spin-dependent scattering on a proton is related to the decomposition 
of the proton spin, in which the greatest uncertainty is the contribution of the strange 
quarks and antiquarks, As''''. Inclusive deep- inelastic scattering experiments favour A^^^ = 
—0.09 ± 0.03 [l9|, whereas analyses of particle production in deep-inelastic scattering are 
quite compatible with As^"* = [20] . In this paper, we compare the neutrino fluxes expected 
for A^P^ = 0, -0.06, -0.09 and -0.12. 

On the other hand, the matrix elements for spin-independent scattering on nuclei are 
related to the vr-nucleon a term, E^jv, whose value is also quite uncertain. The central value 



in a theoretical analysis of low-energy vr-nucleon scattering data is Sttjv = 64 MeV [21 
whereas S^rAr = 36 MeV would correspond to the absence of a strange scalar density in the 
nucleon: (A^|ss|A^) = 0. We consider both these values, as well as the intermediate value 
E^jv = 45 MeV. Lattice calculations are now reaching the stage where they may also provide 



useful information on E^jv [22j , and a recent analysis would suggest a lower value E^jv < 40 
23[. For comparison, we also calculate the flux generated if spin-independent scattering is 
neglected altogether. 

The layout of this paper is as follows. In Section [III we discuss relevant aspects of the 
CMSSM, introducing general aspects of the (mi/2,'n^o) planes for tan/3 = 10 and 55 that 
we use as the basis for our subsequent calculations. In Section [Till we present various astro- 
physical considerations, including a comparison of the LSP capture and annihilation rates 
in the CMSSM, and a discussion of sensitivity to the solar model. In Section [IVl we discuss 
the sensitivities of the annihilation rates to uncertainties in the Standard Model, princi- 
pally the hadronic matrix elements A^'^'* = and Et^jv that control the spin-dependent and 
-independent scattering LSP-matter scattering cross sections, respectively, both of which 
are relevant, and the latter quite important. Then, in Section [Vl we discuss the neutrino 
and neutrino-induced muon fluxes from LSP annihilations, first over the (^«i/2,'"^o) planes 
for tan /3 = 10 and 55, subsequently along the WMAP strips in the coannihilation/rapid- 
annihilation and focus-point strips, and finally for four specific benchmark scenarios. We 
present results for three different neutrino thresholds: 1, 10 and 100 GeV. We find that 
restricted portions of the WMAP strips yields signals that may be detectable in the Ice- 



Cube/DeepCore experiment [2J, |25|. In the coannihilation region, only models with small 
mi/2 are detectable in IceCube/DeepCore. In the case of the focus-point region, a more 
extended part of the strip may be detectable in IceCube/DeepCore for tan/3 = 10. For 
tan /3 = 55, even in the focus-point region only models with small mi/2 may be detectable in 
IceCube/DeepCore. These conclusions are not very sensitive to the uncertainties in the solar 
model and in A^'''* = that controls spin-dependent scattering, but the range of CMSSM 
parameter space that might be observable depends on the hadronic matrix element E^jv 
controlling spin-independent scattering. 



II. CMSSM PARAMETER SPACE 

We set the scene by first discussing the CMSSM parameter space that we explore. Points 
in the CMSSM are in principle characterized by four free parameters, the common gaugino 
mass, mi/2, the common scalar mass, mg, the common trilinear supersymmetry-breaking 
parameter, Aq, and the ratio of Higgs vacuum expectation values, tan/3. There is also 
an ambiguity in the sign of the Higgs mixing parameter, [i: motivated by (7^ — 2 and 
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FIG. 1: The (mi/2,mo) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing regions excluded because there is no consistent electroweak vacuum (dark pink shading), 
or because there is charged dark matter (brown shading), or because of a conflict with 6 — )• 57 
measurements (green shading). Only regions to the right of the black dashed (red dash- dotted) 
line are consistent with the absence at LEP of charginos (Higgs boson). The turquoise strips are 
favoured by the determination of the cold dark matter density by WMAP and other experiments 
fl^ l. and the pink strips are favoured by the BNL measurement of g^ — 2. 



h — )■ 57, we restrict our attention to positive yU ^. The three most relevant parameters for 
this analysis are mi/2,mo and tan/3, and we present our results in the (mi/2,'m.o) planes 
for the trilinear soft supersymmetry-breaking parameter = and two discrete choices 
tan/3 = 10,55. The choice tan/3 = 10 is close to the optimal value we find in a global 
likelihood analysis incorporating all the theoretical, phenomenological, experimental and 
cosmological constraints on the CMSSM parameter space j26[. The choice tan /3 = 55 is 
close the maximum value for which we find consistent solutions of the CMSSM vacuum 
conditions, and represents the most distinct alternative from the favoured case tan /3 = 10. 

In Figure [Hand subsequent figures we display these (mi/2,mo) planes with the following 
constraints implemented: regions lacking a consistent electroweak vacuum are found in the 
upper left corners of the figures (dark pink shading), regions with charged dark matter are 
found in the lower right corners (brown shading), the LEP bound on charginos 27| excludes 
the area to the left of black dashed line and the LEP bound of 114.4 GeV on the Higgs 
mass jl^ is shown by (red) dash-dotted line (lighter Higgs masses occur to the left of this 
line). Here the code FeynHiggs ^] is used for the calculation of m/j. In addition we 
impose agreement with 6 — )■ 37 measurements [30| which excludes the green shaded region 
at relatively low mi/2. We also display the strips favoured by the determination of the 
cold dark matter density by WMAP and other experiments 12| (turquoise shading). Also 



^ We neglect the possibility of significant CP- violating phases in the soft supersymmetry-breaking param- 
eters. 
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shown (shaded pink) is the band favoured by the BNL measurement oi — 2 [3l|], using 
the latest estimate of the Standard Model contribution based on a compilation of e"'"e~ data 
including the most recent BABAR result [32|], which leaves a discrepancy A ((7^ — 2)/2 = 



(24.6 ± 8.0) X 10 |33|] that could be explained by supersymmetry 
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In the case tan/3 = 10 (left panel of Figure [I]), there are two WMAP strips, one close 
to the boundary with charged dark matter where the neutralino coannihilates with the 
lighter stau and other sleptons, and the other close to the boundary of the region with an 
inconsistent electroweak vacuum, the focus-point region. In the case tan /3 = 55 (right panel 
of Figured]), the focus-point strip moves away from the electroweak vacuum boundary at 
large rriQ and small mi/2, and connects with the coannihilation strip. The latter moves away 
from the charged dark matter region at small mo and large mi/2, and bifurcates into strips 
on either side of a funnel where neutralinos annihilate rapidly via direct-channel if, A poles. 

In regions between the WMAP strips and the boundaries set by the absence of charged 
dark matter and the consistency of the electroweak vacuum, the calculated relic density 
falls below the WMAP range [12] • In these regions, for the purposes of our subsequent 
calculations we rescale the halo density by the ratio of the calculated cosmological density 
to the 2-0" lower limit of the WMAP range. The same scaling is applied between the two 
sides of the funnel when tan /3 = 55. Elsewhere, between the WMAP strips, the calculated 
density exceeds the WMAP range if one assumes conventional cosmology, but we neverthe- 
less assume naively that the halo density is saturated by relic neutralinos. This would be 
possible, in principle, if some mechanism increased the entropy density between the epochs 
of dark matter freeze-out, when the temperature T ~ m^/20>5 GeV, and of cosmological 
nucleosythesis, when T < 1 MeV. We do not specifically advocate such a scenario, but dis- 
playing complete (mi/2,'mo) planes enables us to put in context the results we present later 
for models along the WMAP strips. We also present more detailed results for representative 
benchmark points located on these strips. 



III. ASTROPHYSICAL CONSIDERATIONS 

We review here the basic capture and annihilation processes and astrophysical considera- 
tions that affect these processes. The capture and annihilation of LSPs in the Sun depends 
both on the distribution of LSPs in the halo and on the abundances and distribution of 
elements within the Sun. 



For the dark matter halo, we take a non-rotating isothermal sphere [3J] with an rms 
speed of 270 km/s, a disk rotation speed of 220 km/s, and a local dark matter density of 
0.3 GeV/ cm^, and we do not address other halo models in this paper. Our results would scale 
linearly with the local dark matter density as equilibrium between capture and annihilation 
is approached. For our fiducial solar model, we use the AGSS09 model of Serenelli et al. 



[35|, which is based on the proto-solar isotopic abundances found in Ref. [36 . 

Gravitational capture of LSPs occurs when LSPs from the galaxy pass through the Sun, 
scatter off a nucleus, and lose enough energy that they can no longer escape from the Sun's 
gravitational potential. The capture rate per unit stellar mass via scattering off an isotope 



^ The dashed Hnes include the \-a range of /S.{g^ — 2)/2, the sohd hnes the 2-cr range. 
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indexed by i is given by [37|, |38 



^iv) = ^e, fd\l^{u^ + v^)gM, (1) 
aM myjjii J u 

where is the LSP density, m-^ is the LSP mass, rrii is the nuclear mass, o"i is the LSP- 
nuclear scattering cross-section at zero momentum transfer, is the mass fraction of the 
isotope, u is the LSP's velocity far away from the Sun, v{r) is the velocity required for escape 
from the Sun, and 

f{u)^^ldQf{u) (2) 

is the angular-averaged LSP velocity distribution in the galactic neighborhood (in the Sun's 
frame, but outside the solar potential). Here, 



1 pEma-K 

u,v)^-^ diAE)\F,iAE)\'eiP_v'-u'), (3) 



where 

[m^ ± rriiY 

\Fi{/S.E)\^ is a form factor (discussed below), /S.E is the energy imparted to a stationary 
nucleus in a collision with an LSP, and -Emax = ^f^+m^^u^ + f^) is the maximum energy 
that can be imparted by an LSP with velocity Vu"^ + v"^, where y/v? + v^ is the velocity 
of the LSP in the Sun. The quantity Qi{u,v) factors in two effects: (1) only a fraction 
of the scatters result in LSPs losing enough energy to become gravitationally bound, and 
(2) a form factor for finite momentum exchange can suppress higher-energy collisions. In 
the absence of a form factor (l-Fjp = 1), Qi{u,v) is just the fraction of scatters of an LSP 
at velocity y/u'^ + v'^ on a stationary nucleus that yield a gravitationally-bound LSP. The 
theta function 6'(/3_f^ — m^) accounts for the fact that, for sufficiently large LSP velocities, 
scattering into bound orbits is not kinematically possible. The total capture rate is obtained 
by integrating Eqn. ([1]) over the solar profile and summing over the isotopes present in the 
Sun. 

The derivation of Eqn. ([T]) is lengthy and the equations resulting from carrying out the 
integration are long and uniformative, so we do not include them here. Both the source 



of Eqn. ([T]) and integrations are given by Gould in Refs. [37j & j38[. The capture rate is 
arrived at by different methods in those two references; our forms above closely match those 
of the latter, albeit with some small differences ^. We refer the reader to those references 
for technical details and a more thorough discussion of the capture process. Here, we point 
out some of the features of Eqn. ([1]) and Eqn. ([3]): 

• The fraction of scatters that yield capture is largest for nuclei with masses nearest 
that of the LSP. Capture is suppressed for scattering off of light nuclei (mj ^ m^) 
as the LSP loses very little energy in these collisions and is more likely to maintain 
enough speed to escape the solar system and continue back out into the halo. 



The quantity G{u,v) is equivalent to — G{u,v) in Ref. [38|, whilst the LSP-nuclear scattering cross- 
section ai is equivalent to the quantity P+m^rriiaQ Q^. We also differ from Gould in that /(u) is both the 
spatial and velocity phase space density in his derivations (/ d^uf{u) = p^/m^), as opposed to just the 
velocity phase space density here (/ d^w/(u) = 1). 
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• The fraction of scatters that yield capture is higher for scatters that occur deeper in 
the solar potential (larger v), i.e. when closer to the Sun's center. 

• Slow LSPs (small u) require less of an energy loss to become gravitationally bound 
and are more easily captured. 

Scattering in the Sun is complicated by the fact that some LSP-nuclear interactions, 
particularly those involving heavy nuclei, are energetic enough that the finite size of the 
nucleus becomes important. In the zero-momentum-transfer limit, the LSP is expected to 
scatter coherently off the entire nucleus, which can effectively be treated as a point particle; 
an isotropic scattering cross section Uj is assumed in this case. When the inverse of the 
momentum q exchanged in the scatter becomes comparable to or smaller than the size of 
the nucleus, the point-particle approximation is no longer valid and we must include a form 
factor: cTj — )■ (Ji|F(g^)p. We use simple exponential form factors. 



(spin-dependent) 
( spin-indep endent ) 



2 

where is the nuclear mass, is the effective nuclear radius, and AE = ^ is the energy 
lost by the LSP in the collision. Capture rates with a form factor of this form can be found in 



the Appendix of Ref. [37|. We differ here from Gould and other determinations of capture 
rates in that we use a different scale Eq in the exponential for spin-dependent and spin- 
independent scattering For the latter case, this is a reasonable approximation to the more 



accurate Helm form factor for AE ^2Eq |39l. |40| . While the exponential form here becomes a 



poor approximation to the true form factor when AE ^ Eq, the number of scattering events 
at such large energy losses are suppressed to the point that they make little contribution to 
the overall capture rate. The capture rate is thus sensitive mainly to the energy at which 
the form factor becomes important {i.e. Eq), not the shape of the form factor beyond this 
energy. The exponential form factor we use is therefore a reasonable approximation and use 
of e.g. the Helm form factor is not necessary for determining capture rates. This differs from 
direct detection searches, where the signal may be mainly high-energy scatters, for which a 
form factor that is accurate at large AE is necessary. 

As gravitationally-bound LSPs undergo additional scatters on subsequent passes through 
the Sun, they fall inwards and approach thermal equilibrium at the center of the Sun. As 
the population of LSPs increases, the rate at which they annihilate with each other also 
increases; given sufficient time, the annihilation rate will come to equilibrium with the rate 
at which the LSPs are captured. The number of thermalized LSPs follows 

^=C-CeN~ CaN^ , (6) 

where C is the capture rate, Ce parametrizes the evaporation rate of the thermalized LSPs, 
and Ca parametrizes the annihilation rate of the LSPs. Evaporation is only significant for 
LSPs lighter than ~3 GeV and is neglected here, as we do not examine any models with 



^ Note that Gould 37|, |38| defines the exponential scale in terms of the rms radius of the nucleus instead 
of the effective radius used here. 
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such light LSPs. To a good approximation, 

/ 3kT, V'^' 

where < aAV > is the total LSP s-wave annihilation cross section, is the LSP mass, 
and Tc and pc are the temperature and density, respectively, at the center of the Sun. In 
the AGSS09 model, Tc is 1.55 x 10^ K and pc is 151 g/cm^; these two quantities vary by 
negligible amounts for other solar models. Taking Ce = and solving Eqn. ([6]) yields an 
annihilation rate of 

2r = Ctanh2 (^-^^ , (8) 

where is the age of the Sun and r = 1/ ^/CCa- The factor of 2 in 2r is due to the loss of 
two LSPs in each annihilation event. 

It is often assumed that, after the first scatter in which an LSP becomes gravitationally 
bound, it will quickly undergo additional scatters and come to thermal equilibrium in the 
center of the Sun. This assumption is not always valid: in some cases, further scatters may 
occur over time periods longer than the age of the Sun, in which case the LSPs have not yet 
reached thermal equilibrium, or the initial orbits may extend far enough out into the solar 
system that the LSPs are perturbed or gravitationally scattered by the planets, in which 
case they will never come to thermal equilibrium (42], In both cases, the annihilation 



rate will be suppressed. However, for the CMSSM regions examined here, these effects are 
negligible and we ignore them. 



A. Comparison of Capture and Annihilation Rates 



As already mentioned, it might be thought that the LSP density within the Sun is in 
equilibrium, so that the annihilation rate is the same as the capture rate. However, as seen 
in Figure m this is not in general true within the class of CMSSM scenarios discussed here ^. 
In Figure [21 we show contours of the ratio of the annihilation rate to the capture rate. 
Contours where this ratio exceeds 0.8 are in bold for clarity. In particular, for tan /3 = 10, as 
seen in the left panel of Figure [H the LSP annihilation rate is in general much smaller than 
the capture rate, indicating that equilibrium is very poor approximation. The only regions 
in which equilibrium is approached are at very low mi/2 where the LEP Higgs bound (red 
dash-dotted line) is violated, and close to the WMAP strip in the focus-point region, where 
annihilation exceeds 99% of the capture rate. In contrast, the annihilation rate falls to less 
than 1% of the capture rate along the coannihilation strip at large mi/2- The loops centered 
at mi/2 ~ 140 GeV are due to a drop in the relic density due to rapid s-channel annihilation 
through the light Higgs boson. As a consequence, there are drops in both the capture and 
annihilation rates there. The kink and loop in the (^1/2, ""^o) plane for tan/3 = 10 (left panel 
of Figure H]) centered at mi/2 ~ 200 GeV, where the ratio of the annihilation and capture 
rates is relatively high, is due to the enhancement of annihilation just above the W^W~ 
threshold. The fluctuations in the contours of the ratio of the annihilation and capture rates 



^ In calculating the contours in this figure, we use the numerical AGSS09 [35| solar model and the default 
values of the hadronic matrix elements discussed below. 
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FIG. 2: The (mi/2,mo) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing contours of the ratio of solar dark matter annihilation and capture rates, as calculated using 
the AGSS09 model ^shj and default values of the hadronic scattering matrix elements. Equilibrium 
corresponds to a ratio of unity, which is approached for small ?n.i/2 o,nd large tuq. Also shown are 
the theoretical, phenomenological, experimental and cosmological constraints described in the text. 

for mi/2 ~ 400 GeV in the left panel of Figure [2] reflect an increase in the annihilation rate 
above the tt threshold. 

For tan/3 = 55, as shown in the right panel of Figure [2], there are larger regions where 
the annihilation rate exceeds 99% of the capture rate, extending to mi/2 ~ 500 GeV and 
more. In particular, equilibrium is approached in a significant fraction of the coannihilation 
strip, as well as along the focus-point strip, and also in the neighbourhood of the rapid- 
annihilation funnel that extends to mi/2 ~ 1500 GeV, and beyond. This reflects the rapid 
increase in the annihilation cross section as one enters this funnel. We also note that the 
contours of the ratio of the annihilation and capture rates vary rapidly in the neighbourhood 
of the rapid-annihilation funnel, particularly at larger mi/2- As one approaches the funnel, 
the annihilation rate increases, thus increasing the ratio. At higher mi/2, inside the funnel, 
the relic density is suppressed and the rates are scaled. At still larger mi/2 the elastic cross 
section is suppressed and equilibrium is not recovered. 

Being out of equilibrium has important consequences when we examine how variations in 
solar models and, later, scattering cross sections affect the annihilation rates and neutrino 
fluxes. In regions where capture and annihilation are in equilibrium (t© 3> r), 2r oc C so 
that, e.g. a 10% drop in the capture rate would lead to a ~10% drop in the annihilation rate. 
Since, for a given CMSSM model, the neutrino flux is proportional to the annihilation rate, 
the neutrino flux would likewise decrease by ~10%. However, when capture and annihilation 
are well out of equilibrium (t© <C r), 2r oc C^. In this case, a 10% drop in the capture rate 
would lead to a ~19% drop in the annihilation rate and neutrino flux. Thus, the non-linear 
dependence of Eqn. on the capture rate can amplify variations in the annihilation rate 
and hence the neutrino flux relative to the variations in the capture rate found in the various 
cases considered later in the paper. 
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FIG. 3: The {mi/2,mo) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing contours of the ratios of the solar dark matter annihilation rates calculated using the 
GS98 (solid) and A GSSOQph (dot-dashed) solar models to the annihilation rates found using the 
fiducial AGSS09 model pa/ - ^^so shown are the theoretical, phenomenological, experimental and 
cosmological constraints described in the text. 



B. Dependence on the Solar Model 



As discussed above, the rate at which the LSP is captured by the Sun involves determining 
both the rate at which galactic LSPs scatter in the Sun and the fraction of such scatters in 
which the LSP loses enough energy to become gravitationally bound to the Sun. The former 
depends on the LSP-nucleus scattering cross sections, while the latter depends on the mass 
of the nucleus off of which the LSP scatters; both factors may vary significantly among the 
different isotopes present in the Sun. The capture rate is thus sensitive to the abundances 
of elements in the Sun. 

Serenelli et al. 351] have generated models of the Sun for several proto-solar isotopic abun- 
dances^. The models AGSS09 and AGSS09ph are based on the most recent estimates of 
these abundances the latter model differs from the former in that photospheric mea- 
surements instead of meteoritic measurements are used to estimate abundances of refractory 
elements. Solar models based on recent abundance estimates are in conflict with helioseis- 
mological measurements of the Sun; Serenelli et al. provide an additional model (GS98), 
based on older abundance estimates [44], that is compatible with helioseismology. We show 
in Figure [3] comparisons of the annihilation rates in the AGSSOQph and GS98 models with 
the fiducial case (AGSS09) in the planes for tan/3 = 10 and 55. 

The primary difference between the abundance sets and their corresponding solar models 
is an (9(10)% change in the metallicity in the Sun. Whilst these metals represent less 
than 2% of the Sun's mass, we show later that capture via spin-independent interactions 
can be significantly enhanced by heavy nuclei due to a rapid increase in the scattering 



^ We include all 29 isotopes/elements maintained in these models when doing our calculations. 
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cross section. Moreover, the capture of heavy LSPs is favored kinematically by scattering 
on heavier elements. Due to these effects, these trace heavy elements can, in some cases, 
contribute a sizeable portion of the total capture rate. We see in Figure[3]that there is at most 
a 4% difference in the annihilation rates between the models based upon the two most recent 
abundance estimates (AGSS09 and AGSSOQph), so that the discrepancy between meteoritic 
and photospheric measurements is not a significant issue in estimating annihilation rates in 
CMSSM models. However, the model based on older abundance estimates that more closely 
matches helioseismology measurements (GS98) predicts annihilation rates as much as 40% 
higher. In both cases, the discrepancy grows larger as mi/2 increases, since SI scattering, 
which occurs mainly off of the trace heavy elements that differ among the models, becomes 
more important due to the aforementioned effects. In addition, part of the larger variation 
at higher mi/2 can be attributed to non-equilibrium between capture and annihilation as 
discussed previously. For regions noted in Figure [2] as being out of equilibrium, the variation 
in the capture rate for the different solar models is roughly half that of the variation in the 
annihilation rate. For regions in equilibrium, the variation in the capture rate is roughly the 
same as the variation in the annihilation rate. 

It is unclear whether the discrepancies between the recent models and helioseismology is 
due to inaccuracies in the solar modelling or due to the most recent abundance estimates 



being incorrect; see [45[ and references therein for further discussion of this issue. For 
this reason, uncertainties in the solar model still limit the precision to which annihilation 
rates (and, thus, neutrino fluxes) can be predicted for a given CMSSM model, but are not 
important for our conclusions. 



C. Approximating the Solar Potential 

As shown by Gould j38| , the calculation of the capture rate can be simplified if the solar 
potential is approximated by 

*(r) = <fc-^($c-«fs), (9) 

where M(r) is the mass of the Sun interior to a radial distance r, $c and $s are the 
gravitational potentials at the center and surface of the Sun, respectively, and $(r) is the 
potential at an arbitrary point in the Sun. This approximation yields an analytical form for 
the full capture rate, instead of requiring a numerical integration over the solar profile as 
necessary in the exact case. The analytical form is lengthy and can be found in Ref. [ssj. 

The accuracy of the approximation is improved by choosing values for $c and $s such 
that Eqn. ([9]) more closely follows the true potential profile (as determined in the numerical 
models) rather than using the actual central and surface potentials. We use Vc = 1355 km/s 
and Vg = 818 km/s for the central and surface escape velocities (where $(r) = — |t'^(r)), 
which is obtained from a linear fit to the inner 90% (by mass) of the Sun in the AGSS09 
model. The outer 10% is ignored in the fit as the potential can only be well approximated 
by a line for the inner ~90%; inclusion of the outer layer makes the fit much poorer in the 
rest of the Sun. In addition, the dependence of Eqn. ([1]) on the potential generally becomes 
most significant when capture is suppressed at the surface relative to the center; thus, a 
better fit to the inner regions of the Sun is preferred to a fit to the entire solar profile. For 
the parameters given above, Eqn. (Q varies by less than 5% from the actual potential in the 



12 



tan P = 55 , \i>0 




mi/2 (GeV) mi/2 (GeV) 

FIG. 4: The (mi/2,mo) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing contours of the ratio of the dark matter annihilation rates calculated using a full numerical 
integration over the AGSS09 solar profile [35] and using an analytical approximation as proposed by 
Gould. Also shown are the theoretical, phenomenological, experimental and cosmological constraints 
described in the text. 



AGSS09 model for the inner 90% of the Sun. The fits to Vc and Vs for other solar models do 
not change significantly from the above values. 

The Gould approximation requires the further assumption that the abundances of ele- 
ments are uniform throughout the star. This is a reasonable approximation for elements 
heavier than Oxygen, but is not a good approximation for lighter elements. Hydrogen and 
Helium, in particular, are much more and less abundant, respectively, in the center of the 
Sun than elsewhere due to the hydrogen fusion occurring there. Carbon, Nitrogen, and 
Oxygen abundances also vary throughout the star, under the influence of the CNO cycle. 
However, the combined CNO abundance is fairly uniform and, since these three elements 
have similar masses and cross sections, there is little loss in accuracy in taking a uniform 
abundance for them. 

Looking at Figure HJ we see that the numerical integration and the Gould approximation 
agree quite well in general, always within 11% in the tan /3 = 10 plane, and better at high 
mi/2 and in the focus-point region. For tan /3 = 55, the accuracy is even better, with 
the approximation never varying more than 2% from the numerical result. The different 
levels of agreement between the numerical and analytical estimates have two sources: (1) 
varying levels of contribution between the spin-dependent and -independent scattering, and 
(2) amplification of differences when annihation and capture are out of equilibrium. 

The capture rate itself is nearly uniformly overestimated by ~8% for spin-dependent 
interactions and underestimated by ~0.4% for spin- independent interactions when the Gould 
approximation is used. The approximation to the potential, Eqn. OH]), leads to a < 1% 
underestimate. The assumed uniformity of Hydrogen throughout the Sun, however, gives 
a ~8% overestimate of the capture via scattering off of Hydrogen. As Hydrogen accounts 
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for over 99% of the captures via spin-dependent scatters ^, but only a negligible portion of 
the capture via spin-independent scatters, spin-dependent capture is thus overestimated by 
~8%, while spin-independent capture remains underestimated by ~0.4%. Note that Helium 
is not the dominant nuclear target for spin-independent capture, so the poor assumption 
of uniform Helium abundance has little bearing on the results. As seen in Figure |5] and 
discussed in more detail below, spin-dependent scattering is a significant (and dominant 
for tan/3 = 10) contribution to the total capture rate at small mi/2, but is generally less 
significant at higher mi/2, where spin- independent capture takes over. Thus, the Gould 
approximation becomes a better approximation at higher mi/2- 

The total capture rate varies by up to 6% with the approximation, with the largest 
difference occurring for tan/3 = 10, where spin-dependent scattering is more significant. 
Again, non-equilibrium amplifies the variations, increasing 1-6% variations in the capture 
rate to as much as 2-11% variations in the annihilation rate. The largest variation in 
the annihilation rate, which occurs for mi/2 ~ 300 — 500 GeV in the tan/3 = 10 plane 
(left panel of Figure occurs in the region where spin-dependent capture dominates and 
simultaneously annihilation is out of equilibrium. For higher values of mi/2 in the tan /3 = 10 
plane, annihilation is still out of equilibrium, but spin- dependent scattering becomes less 
important. Though the non-equilibrium amplifies variations in the capture rate here, the 
capture rate for spin-independent scattering is much more accurate than the spin- dependent 
case, so the overall annihilation rate in the approximation becomes closer to the true value. 
For the tan/3 = 55 case (right panel of Figure |4]), spin- independent scattering dominates 
everywhere. Because the spin-independent capture rate is very accurately approximated, the 
approximation to the annihilation rate is quite accurate for tan/3 = 55. The shapes of the 
contours in the right panel of Figure H] mainly result from the changing levels of equilibrium 
between annihilation and capture rather than from changes in the relative importance of 
spin-dependent and spin-independent scattering. In the focus-point and funnel regions, 
where annihilation and capture are in equilibrium, the variations in the annihilation rate 
reflect those of the capture rate (~0.5-l%). In the area between the focus point and funnel 
regions, where equilibrium is not maintained, the variation in the annihilation rate grows to 
roughly twice that of the capture rate (~l-2%). 

In principle, numerical integration over the true solar profile is certainly preferable to 
using the Gould approximation, because of its greater accuracy. However, there may be in- 
stances, such as computationally-intensive scans over MSSM parameter space, where avoid- 
ing a numerical integration would be useful ^. For example, using alternative form factors or 
LSP velocity distributions may require the integrals in Eqn. ([1]) to be evaluated numerically, 
which could make determining the capture rate a computationally-slow process where reduc- 
ing the number of numerical integrations would be beneficial time-wise. Also, we note that 
using this approximation introduces inaccuracies no worse than those due to choices in solar 
models and uncertainties in hadronic parameters, as discussed below. There is also room 
for improvement in this approximation as the chief source of its inaccuracy — the assumed 
uniformity of the Hydrogen abundance — could be significantly alleviated by, e.g., breaking 



We have included spin-dependent scattering on all significant isotopes in the Sun (notably Helium-3 and 
Nitrogen- 14), but find that elements other than Hydrogen account for much less than 1% of the capture 
from this type of scattering in the models considered here. 
® The same could be accomplished in some cases — but not all — by tabulating the integration, as done in 
DarkSUSY j46|. 
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-0.09 ± 0.03 


[19] 



TABLE I: Strong-interaction parameters used to determine neutralino-nucleon scattering cross sec- 
tions, with estimates of their experimental uncertainties. As discussed in the text, the most impor- 
tant uncertainties are those in cq, and /\'f\ We use as defaults for them the central values 
given in the Table. 

Eqn. ([9]) into two parts, corresponding to the inner and outer portions of the Sun, and using 
more appropriate Hydrogen abundances in each part. 

IV. STANDARD MODEL PARTICLE PHYSICS UNCERTAINTIES 

We now turn to the uncertainties in the annihilation rates induced by the uncertainties in 
the hadronic matrix elements that enter into dark matter scattering rates, which are listed in 
Table [H They include those in the quark masses, expressed as md,c,b,t and the ratios m^jmA 
and m^jm^^ those in the matrix elements l\N\qq\N) - which are related to the change in the 
nucleon mass due to non-zero quark masses, denoted by ctq, and therefore to the vr-nucleon 
cr term S^^ as discussed later - and the axial-current matrix elements (iV| 57^755 |A^), which 

are related to the quantities A^^^ ag'''' and also discussed later. The uncertainties 

in the elastic scattering cross section induced by the uncertainties in the quark masses, 
apart from the top quark, are negligible. However, cross-section uncertainties induced by 
the uncertainties in the matrix elements (A^|^7^75q'|iV) and (iV|^q'|A^) are important, as 
we discuss below. In particular, the uncertainties induced by our ignorance of A^P^ and 
particularly E^^r should not be neglected. 

As already mentioned, it might be thought that the dominant dark-matter scattering 
process inside the Sun is spin-dependent scattering on the protons that comprise most of the 
solar mass. However, this is not necessarily the case in the CMSSM, as seen in Figure [5] which 
shows the annihilation rate when only spin-dependent scattering is included, relative to the 
total annihilation rate when both spin- dependent and -independent scattering is included. 
We see in the left panel for tan/3 = 10 that, whereas spin-dependent scattering is indeed 
dominant at small mi/2 and large mg, this is a poor approximation already in the portion of 
the coannihilation strip that is favoured by (7^ — 2, and spin-independent scattering actually 
dominates at large mi/2 along the coannihilation strip. For tan/3 = 55 (right panel), spin- 
independent scattering dominates across the plane and signicantly so in the coannihilation 
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FIG. 5: The (mi/2,mo) planes for the CMSSM for Aq = and (left) taii/3 = 10, (right) tan/3 = 55, 
showing contours of the ratio of the solar dark matter annihilation rate calculated using only spin- 
dependent scattering to the total annihilation rate including also spin-independent scattering. Also 
shown are the theoretical, phenomenological, experimental and cosmological constraints described 
in the text. 

strip and in the rapid- annihilation funnel ^ . The increase in the relative importance of 
spin-independent scattering at large mi/2 may be traced to the importance of the reduced 
neutralino mass in the kinematics of the scattering process. As the neutralino mass 
increases, the dependence on disfavours capture via scattering off light nuclei such as 
the proton that dominates spin-dependent scattering, and increases the relative importance 
of spin- independent scattering, which is dominated by scattering off more massive nuclei. 
Furthermore, because the spin-independent cross section increases with tan /3, its importance 
is amplified in the right panel for tan/3 = 55. As in previous figures, we see in the right 
panel of Figure for tan /3 = 55 that the contours of the ratio are relatively close together 
in the neighbourhood of the rapid-annihilation funnel. In this region the annihilation rate 
varies rapidly and approaches equilibrium with the capture rate. Without spin-independent 
scattering, the capture rate there would be reduced, driving annihilation out of equilibrium 
and suppressing the annihilation rate and the resulting neutrino flux. 

The results shown in Figure [5] imply that one must discuss the uncertainties in both spin- 
dependent and -independent LSP scattering, as we discuss in more detail in the subsequent 
subsections. 



^ In calculating the contours in this figure, we use the default values of As and S^jv discussed below. 
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A. Sensitivity to Spin-Dependent Scattering 

The spin-dependent (SD) part of the elastic x-nucleus cross section can be written as 

asD = —GlmlA^JiJ+l), (10) 

TT 

where rrir is the reduced neutrahno mass, J is the spin of the nucleus, 

A = i(ap(Sp) + an(S,)) , (11) 



{Sp) and {Sp) are the spin content of the proton and neutron groups, respectively (53|, and 

a(p) , _ "2g ^ (n) 



«P = E ^A?^ «n = E . (12) 



The factors Ag^"* parametrize the quark spin content of the nucleon and are significant only 
for the light (u,d,s) quarks. A combination of experimental and theoretical results tightly 
constrain the linear combinations 148| 



^(P) ^ ^(p) _ ^(P) _ ^ 2694 ± 0.0028 (13) 

and |5ll|52 

^(P) ^ ^ip) ^ ^(P) _ 2^{P) = 0.585 ± 0.025. (14) 

However, determination of the individual A[^^ requires a third piece of information, usually 
taken to be the strange spin contribution As'''', as extracted from inclusive deep- inelastic 



lepton-nucleon scattering. Using, e.g., the recent COMPASS result [19|, one has 



A(p) = -0.09 ± 0.01 (stat.) ± 0.02 (syst.) ^ -0.09 ± 0.03 , (15) 

where, conservatively, we have combined linearly the statistical and systematic uncertainties. 
Using this range for As'''' , we may express A^^'' as follows in terms of known quantities: 

A(P) = ^ (^a'f^ + ) + A(P) = 0.84 ± 0.03 , (16) 
A^P) = ^ (4") - a^P)) + A(P) = -0.43 ± 0.03 . (17) 

The above two uncertainties are almost completely correlated with that of As'''', however, 
and the uncertainties in 03''^ and a^g^ are negligible by comparison. We use the central value 
in Eqn. (fT5|) as our default, namely Ai^^ = —0.09 

For comparison, in Figure [6] we compare the annihilation rates calculated using this 
default value and with A^'''' = —0.12 and —0.06. We recall, however, that analyses of vr and 
K meson production in deep-inelastic scattering have been interpreted as suggesting that 



10 



We recall that the proton and neutron scalar matrix elements are related by an interchange of Au and 
Ad, z.e., a(") = A^P\ a(") = aIF'> and A^'^) = A^^^ 
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mi/2 (GeV) mi/2 (GeV) 



FIG. 6: The (w.1/21 ''tt-o) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing contours of the ratios of the annihilation rates calculated assuming A^P^ = 0.00 (dotted 
lines), —0.06 (solid lines) and —0.12 (dashed lines) to calculations with the default value = 
—0.09. Also shown are the theoretical, phenomenological, experimental and cosmological constraints 
described in the text. 



As may be compatible with zero [20] . Accordingly, we also present in Figure |6] results for 
this extreme value of As''''. 

In general, the spin- dependent scattering rate increases and decreases with A^'''*, and 
the effect is quite symmetric for As'''' G (—0.12,-0.06). However, the effect is quite small 
for A^''^ within this range, differing from the rate calculated with our default assumption 
As''^ = —0.09 by at most ~ 10% for small mi/2 when tan/3 = 10 (left panel of Figure E]) and 
considerably less for tan /3 = 55 (right panel of Figure [6]) However, the reduction may be 
considerably larger if A^'''' = 0, potentially exceeding 20% at small mi/2 when tan/3 = 10. 



B. Sensitivity to Spin-Independent Scattering 

The spin-independent (SI) part of the cross section for scattering on a nucleus {Z, A) can 
be written as 

asi = ^[Z/, + (A-Z)/„f , (18) 
7r 



where is the x-nuclear reduced mass and 

fN _ AN) «3g 2 (TV) a?,q 



rriN ^ m„ 27"" ^ m„ 

5=u,d,s ^ g=c,b,t ^ 



The contours near the rapid-annihilation funnel again reflect the equilibrium/non-equilibrium effects on 
the annihilation rate discussed earlier. 
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for = p or n. The parameters Z^^'' are defined by 

m^4f = {N\m,qq\N) = m,Bf\ 



where 154 55 



AN) 

JTG 



(N) 



1- 

5=u,d,s 

We take the ratios of the light quark masses from ^ 



TTi TTl 

— = 0.553 ±0.043, — = 18.9 ±0.8, 



(20) 
(21) 

(22) 



and the other quark masses are taken from j48|, except for the top mass, which is taken from 
the combined CDF and DO result These masses, as well as other experimental quantities 
that will arise in the calculation of the hadronic matrix elements, appear in Table [B 
Following |56|, we introduce the quantity: 



r(p) r(p) 



B. 



(p) 



(P) 



1.49, 



(23) 



which has an experimental error that is negligible compared with others discussed below, 
and the strange scalar density 



2B, 



(JV) 



y 



i?f)±i?r' 



(24) 



In terms of these quantities, one may write 



B 



(p) 



2+{{z-l)xy) 



(25) 



^(P) 2x z-{{z-l)xy) ' 
Proton and neutron scalar matrix elements are related by an interchange of i?u and B^, i.e., 



B^^)=b[?\ # = R(p) and Bi'^'^ = Bi^\ 



The TT-nucleon sigma term, E^jv, niay be written as 

1 



2("^u±md)x (Sf)±sf)) , 



and the coefficients /t, may be written in the forms: 

2S. 



rriN 
rriN 



TTtu 

2S 



-^(i + t;)(i + f)' 



^^(i + tf)(i + f) ' 



l^TTJv y 



rnu\ 1 
"Id ^ 



(26) 
(27) 

(28) 
(29) 

(30) 
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FIG. 7: The {'mi/2,'mo) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan^ = 
55, showing contours of the ratios of the solar dark matter annihilation rates calculated using 
SttAt = 36 MeV (dashed lines) and 45 MeV (solid lines) to calculations with the default value 
SttAt = 64 MeV. Also shown are the theoretical, phenomenological, experimental and cosmological 
constraints described in the text. 



where we have dropped the (A^) superscript from /t^ and B, 

The effects of the uncertainties in the /t, were considered in [57[, and we were motivated 
to study jl6[ variations in the value of y by re-evaluations of the 7r-nucleon sigma term T^-^n, 
which is related to the strange scalar density in the nucleon by 

y=l- o-q/S^jv • (31) 
The value for ctq given in Table |T]is estimated on the basis of octet baryon mass differences to 



be (To = 36 ±7 MeV |50l. |58H60| . Based on recent determinations of S,rjv at the Cheng-Dashen 
point [2l|, one finds 

E^jv = (64 ± 8) MeV . (32) 

We use the central value in Eqn. f p2|) as our default, namely St^jv = 64 MeV. 

In Figure [71 we compare the annihilation rates calculated with this default value and with 
S^jv = 45 or 36 MeV, the latter corresponding via Eqn. fl3T]) to the central value for ctq, i.e., 
y = 0. As could be expected from the previous analysis demonstrating the importance of 
spin-independent scattering at large mi/2, we see that the annihilation rate decreases very 
substantially if S^jv is decreased to 45 or 36 MeV As in the case of direct dark matter 
detection via spin-independent scattering, it is also very important for the interpretation of 
indirect dark matter searches via neutrinos from dark matter annihilation inside the Sun to 
pin down the magnitude ofJ^j^f^. 



The contours near the rapid-annihilation funnel again reflect the non-equilibrium effects on the annihila- 
tion rate discussed earlier. 
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V. NEUTRINO AND MUON FLUXES 



Neutrino detectors such as IceCube/DeepCore [2J, |25[ aim to detect the capture of LSPs 
in the Sun using neutrinos produced in the LSP annihilations. The primary experimental 
signature is the passage of neutrino-induced muons through the detector. Here, we examine 
the flux of neutrinos and neutrino-induced muons in the tan /3 = 10 and 55 planes, as well 
as along the WMAP-preferred focus-point, coannihilation, and funnel-region strips (the last 
for tan/3 = 55 only). We also study the neutrino and flux spectra for several benchmark 
scenarios given in Table HTl 

Determining the neutrino spectrum in a detector involves two steps: (1) the production 
of neutrinos from LSP annihilations, and (2) propagation of the neutrinos from the interior 
of the Sun to the detector. Low energy neutralinos do not generally annihilate directly into 
neutrinos; instead, neutrinos are produced in decays/showers of the primary annihilation 
particles such as W bosons, top quarks, or tau leptons {e.g. XX ~^ 't^i with r — )• /iz/^z/,-). 
Determination of the neutrino spectra in such showers is complicated by the fact that they 
occur in the dense center of the Sun, so that the primary particles may lose energy before 



decaying |6 11464 



Once neutrinos are produced, they must travel through the Sun, where they may un- 
dergo charged-current or neutral-current interactions that absorb the neutrinos or reduce 
the neutrino energies, respectively. The Sun is mainly transparent to neutrinos well below 
~100 GeV, but becomes opaque when the energies reach ~200-300 GeV; thus, high-energy 
neutrinos are heavily suppressed. Oscillations between neutrino species between the Sun 



and Earth must also be taken into account 65 



We use the neutrino and neutrino-induced muon spectra determined by WimpSim [6 



(and used within DarkSUSY 46|), which simulated both of the steps described above to gen- 
erate spectra for several LSP annihilation channels. Neutrino production and propagation 
has also been simulated in Ref. j63| 

For our neutrino spectrum results here, we neglect Higgs annihilation channels and include 
only annihilations into quarks, leptons, W and Z bosons. The Higgs channels account for less 
than 3% of the annihilations in all the tau/S = 10 region we have considered and, apart from 
a small region at mi/2 > 1900 GeV just above the charged dark matter constraint (brown 
region in the planar figures), Higgs channels account for less than 10% of the annihilations in 
the tan /3 = 55 plane. In most cases, and all of the WMAP strips, the Higgs channels account 
for ~1% or less of the annihilations. The error introduced by neglecting the Higgs channels 
is thus far smaller than variations induced by the choice of solar model or uncertainties in 
the hadronic parameters discussed previously. 

The muon energy detection threshold i?th is an important consideration in determining 
the sensitivity of neutrino detectors to annihilations in the Sun. Increasing the exposure area 
of a detector often requires a sacrifice in the energy threshold. Since high-energy neutrinos 
from the Sun (and thus muons) are suppressed, high thresholds may inhibit detectors from 
observing this type of signal. The full IceCube detector will have a fairly large effective area, 
but gradually loses sensitivity to muons below energies ~100 GeV due to the large spacing 
between strings [2^]. A smaller portion of the detector, referred to as DeepCore, will be 
instrumented more densely and with higher efficiency phototubes to allow for detection of 



During the preparation of this paper, an error was discovered in the anti-neutrino fluxes of Ref. [67.] that 
prevented us from making fuh use of their resuhs at this time. 
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lower-energy muons; DeepCore can potentially detect muons down to energies of 10 GeV 
or lower j25|. The effective neutrino detection area of the combined IceCube/DeepCore 
detector falls from ~1 m^ at 1 TeV, to ^10'^ m^ at 100 GeV, to ^10"^ m^ at 10 GeV, with 
IceCube providing nearly all of the effective area above 100 GeV and DeepCore providing 
nearly all of it below 30 GeV. 

The ability of IceCube/DeepCore to observe a signal above background for a particular 
CMSSM model depends on the neutrino spectrum, angular spread of the induced muons, 
backgrounds, and detector geometry 68, 6^. The predominant backgrounds are muons 



and muon neutrinos generated from cosmic ray interactions in the atmosphere; cosmic ray 
interactions in the Sun may also contribute to the background [70|. When looking for a 
signal from a point source such as the Sun, the diffuse background is reduced by making 
angular cuts on the observed muon event directions. A signal analysis is complicated by 
the fact that the angular spread of solar neutrino-induced muons is dependent on neutrino 
energies, so that optimal data cuts are dependent on the expected neutrino spectrum for a 
particular CMSSM model. In addition, the angular and energy resolutions and the effective 
area for muon events, particularly at low energies, is dependent on the detector geometry. 
Such a careful analysis is beyond the scope of this paper. 

Still, we can roughly approximate the IceCube/DeepCore sensitivity to estimate the de- 
tectability of various CMSSM models. IceCube/DeepCore has generated a conservative 
sensitivity limit for muons above 1 GeV in their detector assuming a particular hard anni- 
hilation spectrum [g^ (rf for < 80 GeV, W^W~ at higher masses). The fraction of 
annihilations into channels that yield hard spectra {xx ~^ ZZ,W~^W~ ,tt or rr) is never 
trivial in the CMSSM models we consider in this paper. In the tan /3 = 10 case, 70-99% of 
the annihilations are into one of these hard channels, whilst in the tan /3 = 55 case 20-50% 
of the annihilations are into one of these channels. The fractions of annihilations into soft 
spectra channels {e.g. xX ~^ ^) ^-^e 1-30% and 50-80% for the tan/3 = 10 and 55 cases, 
respectively. Thus, in all cases, the spectrum produced in a CMSSM model will contain 
a significant hard component. Since these CMSSM neutrino spectra differ from the spec- 
tra used by IceCube/DeepCore to generate their muon flux limits, the IceCube/DeepCore 
muon flux limits cannot be directly applied to particular CMSSM models. However, the 
significant hard component of the CMSSM spectrum should be similar enough to the hard 
spectrum channels analyzed by IceCube/DeepCore that the sensitivity limits are within a 
factor of 2 (tan/3 = 10) or 2-5 (tan/3 = 55) of the sensitivity that would be obtained had 
the actual neutrino spectra been used in the analysis. Hence, in the tan/? = 55 case, the 
limits obtainable are expected to be somewhat weaker than in the tan /3 = 10 case, due to 
the lower fraction of annihilations into hard channels. 

The IceCube/DeepCore projected sensitivity is shown in Figures [9lfTT| normalized to the 
total expected muon flux above 1 GeV (though the experiment is only weakly sensitive to 
muons at such low energies) We note that the IceCube/DeepCore collaboration charac- 
terizes this as a "conservative" limit, and that the sensitivity could be further improved in 
several possible ways. For example, (1) data are generally used only from periods when the 
Sun is below the horizon, as this reduces the cosmic-ray backgrounds. However, it may be 
possible to use the surrounding IceCube portion of the detector to veto cosmic-ray- induced 
events in the DeepCore portion of the detector, allowing data from daylight hours to be 



If the actual spectrum were to be entirely soft, the IceCube/DeepCore muon flux sensitivity would be 
weakened by an order of magnitude or more. 
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included in the analysis and hence potentially doubling the exposure at low energies. In 
addition, (2) even a few more strings placed in the DeepCore region (a very interesting pos- 
sibility in light of our analysis) would significantly improve sensitivity to low-energy muons. 
For these reasons, the IceCube/DeepCore muon flux sensitivity limits should be taken only 
as an order of magnitude estimate when assessing below the detectability of the CMSSM 
models below. 

A. Fluxes in (m,i/2,?7io) Planes 

We present in Figure E] contours of the neutrino-induced muon fluxes in the (mi/2,wio) 
planes for tan /3 = 10 (left) and tan /3 = 55 (right), for three values of the detector threshold: 
1 GeV (top), 10 GeV (middle), and 100 GeV (bottom). The other phenomenological, exper- 
imental and cosmological constraints are as in Figure |2j As the aforementioned sensitivity 
limits cannot easily be applied to these figures, we can expect, to a very rough approxima- 
tion, IceCube to detect muon fluxes on the order of 10 or 10^ /km^/yr above ~100 GeV and 
DeepCore to detect muon fluxes on the order of 10^ or 10'^ /km^/yr above ~10 GeV. 

For tan /3 = 10, the muon flux attains its maximum value close to the focus-point region. 
The picture does not change much when the detector threshold increases from 1 to 10 GeV, 
but for threshold 100 GeV the contours show lower muon fluxes. IceCube/DeepCore is in 
this case expected to probe all the focus-point region and much of the bulk region. Note 
that although the highest flux contour we display is 500 /km^/yr, the flux continues to 
increase as one approaches the region where radiative electroweak symmetry breaking no 
longer is viable. The focus-point strip is very close to this boundary and indeed the flux does 
exceed 10^ /km^/yr when the relic density is in the WMAP range. This will be seen more 
clearly when we discuss the the WMAP strips in the next subsection. The focus-point region 
produces the largest fluxes observable by IceCube alone (high threshold, bottom left panel 
of Figure [H]); these high-energy muon fluxes may even be detectable without including the 
sensitivity to lower energy muons that DeepCore provides. The same holds for tan /3 = 55: 
the neutrino-induced muon flux is again largest in the focus-point region, close to the region 
theoretically excluded due to the absence of electroweak symmetry breaking. However, the 
fluxes are likely too low to be observed by IceCube/DeepCore, except for the lowest values 
of mi/2 along this strip. The cosmologically-favoured region around the rapid-annihilation 
funnel yields much smaller values for the muon fluxes, below the expected experimental 
sensitivities. 

B. Fluxes along WMAP Strips 

In Figure HI we plot the neutrino-induced muon fluxes in the cosmologically-favored 
regions of the CMSSM that are also compatible with the rest of the phenomenological 
and experimental constraints, namely the regions that we call WMAP strips. These strips 
are determined by fixing the value of ttiq so that at each value of mi (or equivalently 
m-^ ~ 0.43mi/2) the relic density is within the 2 a WMAP range [l2|. We pick again 
the same representative values of tan/3 = 10 (left) and 55 (right). For tan/3 = 10, there 
are two such strips corresponding to either the co- annihilation or focus-point regions. For 
tan/3 = 55, the co-annihilation strip runs into the funnel region at large m^. At very large 
m^, we have sampled both sides of the funnel resulting in a 'doubling back' of the curves 
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FIG. 8: The (1^.1/2,1^0) planes for the CMSSM for Aq = and (left) tan/3 = 10, (right) tan/3 = 55, 
showing contours of the neutrino-induced muon flux through a detector with thresholds of 1 GeV 
(top), 10 GeV (middle), and 100 GeV (bottom). Also shown are the theoretical, phenomenological, 
experimental and cosmological constraints described in the text. 
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FIG. 9: The CMSSM muon fluxes though a detector calculated for Aq = Q and (left) tan /3 = 10, 
(right) tan/3 = 55, along the WMAP strips in the coannihilation/ funnel regions (solid) and the 
focus-point region (dashed). Fluxes are shown for muon energy thresholds of (top to bottom) 1 GeV, 
10 GeV, and 100 GeV. Also shown is a conservative estimate of sensitivity of the IceCuhe/ Deep Core 
detector (dotted), normalized to a muon threshold of 1 GeV, for a particular hard annihilation 
spectrum (tt for < 80 GeV, W~^W~ at higher masses). The IceGuhe/DeepCore sensitivity 
shown does not directly apply to the CMSSM flux curves, hut can he treated as a rough approximation 
to which CMSSM models might be detectable, as discussed in the text. 

we display. 

A conservative IceCube/DeepCore sensitivity limit for muon fluxes above 1 GeV is also 
shown in the figure. As discussed previously, this sensitivity cannot be directly applied 
to the CMSSM fluxes shown as the spectra differ. However, the limit should be within a 
factor of 2 or so of the muon fluxes necessary to detect these CMSSM models (recall that, 
if additional strings are added to DeepCore and/or data can be used from periods when the 
Sun is above the horizon, the IceCube/DeepCore sensitivity could be as much as an order 
of magnitude better than shown). 

As observed in Figure [HI the muon flux is largest in the focus-point region for both values 
of tan/3. In Figure [9] (left) we plot the values for the muon fluxes along the coannihilation 
strip (solid curve) and the focus-point strip (dashed curve) for tan/3 = 10. The blue, red 
and green curves correspond to threshold energies 1, 10 and 100 GeV. In the focus-point 
region, the muon flux is ~ 10^ - 10'^ /km^/yr, making detection in IceCube/DeepCore a 
likely possibility, particularly for lighter LSPs (smaller mi/2)- For tan/3 = 55 (right plot), 
the solid curve represents both the coannihilation and funnel strips. The double covering 
of ~ 550 — 750 GeV is due to the shape of the WMAP strip in the funnel region, 
as noted above. As already observed, the values of the muon fluxes are similar in the 
coannihilation/funnel region and in the focus-point region for tan/3 = 55, but are smaller 
than those for tan/3 = 10, especially in the focus-point region. 

We see in Figure [TO] that the muon fluxes are quite insensitive to the value of A^'''* within 
the (plausible) range studied. On the other hand. Figure [11] shows that the muon fluxes are 
more sensitive to the value of St^jv, particularly for tan/3 = 55 and along the coannihilation 
strip for tan/3 = 10 at larger mi/2- The comparisons shown are for a threshold energy of 
1 GeV, but similar conclusions apply for thresholds of 10 and 100 GeV (not shown). The 
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FIG. 10: The CMSSM muon fluxes though a detector calculated for Aq = and (left) tan/3 = 10, 
(right) tan/3 = 55, along the WMAP strips in the coannihilation/ funnel regions (solid) and the 
focus-point region (dashed). The fluxes are shown for A^^^ values of 0, —0.06, —0.09, and —0.12 
with a detector threshold of 1 GeV. Also shown is a conservative approximation of sensitivity of 
the IceCuhe/ Deep Core detector (dotted), as described in FigurelM o-nd the text. 
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FIG. 11: The CMSSM muon fluxes though a detector calculated for j4o = and (left) tan/3 = 
10, (right) tan/3 = 55, along the WMAP strips in the coannihilation/ funnel regions (solid) and 
the focus-point region (dashed). The fluxes are shown for S^^jv values of (top to bottom) 64, 45, 
and 36 MeV with a detector threshold of 1 GeV. Also shown is a conservative approximation of 
sensitivity of the IceCube/DeepCore detector (dotted), as described in Figure\^and the text. 



value of S^jv clearly impacts the ability of IceCube/DeepCore to detect tan/3 = 55 models. 



C. Fluxes for Benchmark Points 

We now display complete spectra for the CMSSM benchmark scenarios C, E, L and M 
of [isj. These are representative, respectively, of the low-mass region of the coannihilation 
strip, the focus-point region, the coannihilation region at larger mi/2 and tan/3, and the 
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TABLE II: CMSSM parameters and results for benchmark models C, E, L, and M of JlA l. Eth 
is the neutrino /muon energy threshold. Rates and fluxes are determined using the AGSS09 solar 
model and the central values found in Tahle\^for the hadronic parameters. 

rapid- annihilation funnel region. The benchmark parameters, capture/annihilation rates, 
and fluxes are given in Table HTl Models E, L, and M are in or nearly in equilibrium between 
capture and annihilation, whereas model C is somewhat out of equilibrium. 

The differential neutrino fluxes for these scenarios are shown in Figure [121 The highest 
differential flux is found in model E, partly because it has a larger annihilation rate than 
the other models and partly because its neutrino spectrum is spread over a smaller range 
of energies (at 117 GeV, this model has the smallest LSP mass of the four). Model M has 
the lowest differential flux, but the spectrum is spread out over a wider range of energies; 
indeed, the spectrum for this model continues beyond the region shown. In all four models, 
the spectra end at the corresponding LSP mass. The plateau seen in model E arises from 
W — )■ /iz/^ decays in the dominant XX ~^ WW annihilation channel. Since the LSP mass is 
not much higher than the W mass, the Ws are produced with relatively low speeds. In this 
case, the neutrinos from the W decays are kinematically limited to fall within the energy 
range corresponding to the plateau in the spectrum. 

Figure [12] shows that the neutrino flux (solid) is slightly smaller than the anti-neutrino flux 
(dashed), particularly at higher energies. This arises from the higher neutral- and charged- 
current scattering cross sections for neutrinos, compared to anti-neutrinos, when passing 
through the Sun. The difference between neutrino and anti-neutrino fluxes grows more 
pronounced at higher energies. However, both neutrinos and anti-neutrinos are suppressed at 
higher energies, due to the increasing opacity of the Sun: neutrinos are just more suppressed 
than anti-neutrinos. 

Figure [13] shows the differential muon fluxes from these neutrino spectra. Model E still 
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FIG. 12: The neutrino spectra in a terrestrial detector for (top to bottom) benchmark scenarios E, 
L, C, and M fl^J. Neutrinos (solid) and anti-neutrinos (dashed) are shown separately. 
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FIG. 13: The muon spectra in a terrestrial detector for (top to bottom) benchmark scenarios E, L, 
C, and M Muons (solid) and anti-muons (dashed) are shown separately. 

has the highest differential flux, but model M has increased relative to the other models over 
the intermediate energies shown in the figure, as compared to Figure [121 This is because the 
higher-energy neutrinos in this model are capable of inducing additional muons at a range 
of lower energies. 

Whilst model E has the largest muon flux, it is all at relatively low energies that IceCube 
alone cannot detect. With a flux of just under 10'^ /km^/yr above 10 GeV (see Table HTll . 
however, its flux lies above the expected DeepCore detection level. Model L has the next 
highest spectrum, but would generate ~50 muons/km^/yr above 10 GeV, probably falling 
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just below the DeepCore detection level. However, its spectrum extends up to 190 GeV, 
yielding a significant flux above 100 GeV of 2.7 /km^/yr. While not high enough to be 
detected by IceCube alone, it is not far below the expected IceCube detection level either. 
The combination of DeepCore's low energy sensitivity and IceCube's large volume make 
IceCube/DeepCore perhaps able to detect this model. Model M, while producing high- 
energy muons to which IceCube is most sensitive, has fluxes that are too low to be detected. 
Model C has neither the total flux nor an energy spectrum favorable for detection. 



VI. CONCLUSIONS 

We have studied in detail the observability of muons produced by neutrinos and antineu- 
trinos generated by the annihilations of CMSSM LSPs trapped inside the Sun, focusing 
in particular on models lying along the parameter strips where the relic LSP density falls 
within the range favoured by WMAP [12] . 

For various reasons, discussed in Section |Vl it is not straightforward to apply direct de- 
tectability limits for IceCube/DeepCore experiments on the CMSSM parameter space. The 
main complication is that the ability of the experiment to observe a signal above background 
depends on the shape of the neutrino spectrum. Nevertheless, a conservative estimation of 
the IceCube/DeepCore detectability limits can be made, that being 10^ to 10^ /km^/yr for 
the muon flux above 1 GeV. Based on this, our conclusions are not very encouraging: along 
the coannihilation strips, we find that only models near the low-mi/2 ends of the strips for 
tan /3 = 10 and 55 may be detectable in IceCube/DeepCore. Nearly all of the focus-point 
strip may be detectable by IceCube/DeepCore for tan/3 = 10; some of this strip may ac- 
tually be detectable with the bulk {i.e. non-DeepCore) IceCube detector region alone. For 
tan /3 = 55, IceCube/DeepCore will only be sensitive to the low-mi/2 end of the focus-point 
region. 

As for the benchmark scenarios, we find that model E (from the focus-point region) 
is expected to be detectable in IceCube/DeepCore, mainly due to the low energy muon 
sensitivity of the DeepCore component as this model produces few high energy muons that 
IceCube alone is sensitive to. On the other hand, the muon spectrum in model L (in the 
coannhilation region) extends to an energy of 190 GeV, and both the DeepCore and bulk 
IceCube components contribute to making this model borderline detectable. Model M (in the 
rapid-annihilation region), while also producing high-energy muons, produces too few to be 
detectable in IceCube/DeepCore, and model C (near the low-mi/2 end of the coannihilation 
strip) has neither the flux nor the energy spectrum suitable for detection. These conclusions 
are quite robust with respect to uncertainties in the solar model. However, we do note 
that the calculated muon fluxes are quite sensitive to the assumed value of the 7r-nucleon a 
term St^jv, which controls the magnitude of spin- independent LSP scattering: we repeat our 



previous plea [17[ to our experimental colleagues to reduce the uncertainty in St^jv- 

These studies indicate that the indirect search for dark matter via annihilations that yield 
high-energy neutrinos and hence muons may not be the most promising way to discover 
supersymmetry, at least within the restrictive CMSSM framework. On the other hand, 
there are representative models where the calculated fluxes are close to the sensitivities 
of DeepCore and/or IceCube. Therefore, if other evidence for superymmetry were to be 
found, e.g.., either at the LHC or in the direct search for dark matter scattering, it would 
be interesting to increase the sensitivities of DeepCore and IceCube, e.g., by decreasing the 
threshold or by increasing the effective surface area. Une affaire a suivre! 
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